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■ Abstract 

We point out that the temperature fluctuations of cosmic microwave background 
(CMB) can be generated in a way that is different from the one usually assumed in 
slow-roll inflation. Our mechanism is based on vacuum fluctuations of flelds which 
are at rest at the bottom of the potential, such as Kaluza-Klein modes or string 
^ I excited states. When there are a large number (typically of order N ~ 10^^) of fields 

■ with small mass in unit of Hubble parameter during the inflationary era, this effect 
can give significant contributions to the CMB temperature fiuctuations. This number 
N makes it possible to enhance scalar perturbation relative to tensor perturbation. 
Comparison with the observed amplitudes suggests that models with string scale of 
order 10~^ of 4D Planck scale are favorable. 
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1 Introduction 



Observation of cosmic microwave background (CMB) provides an excellent opportunity 
for testing theories of high energy physics. The CMB radiations are the photons emitted 
at the era of recombination reaching us almost unscattered. It has a very homogeneous 
distribution over the whole sky with thermal spectrum at T ~ 2.7K with fluctuations 5T /T 
of order 10~^. Temperature fluctuation is directly related to the gravitational potential 
at the last scattering surface by the relation, 5T/T = — $/3 (see e.g. [1]). Gravitational 
potential $ is essentially frozen in the matter or radiation dominated universe, thus the 
observation of CMB enables us to trace back the universe to the era much earlier than 
recombination. 

It is believed that there has been a period of exponential expansion (inflation) in the 
early universe [2]. Had the universe been decelerating (matter or radiation dominated) 
since the beginning, the observable universe would have to be made of many spatial regions 
which have been initially independent, making it difficult to explain the homogeneity of 
our universe. Exponential expansion brings these regions in causal contact in the past. 
This is the only compelling resolution of this horizon problem. 

The fluctuations generated during inflation has nearly scale invariant spectrum. At 
each time 6t ~ (where H is the Hubble parameter of inflation), fluctuations of order 
6(j) ^ H will be created in a spatial region of horizon size ~ H~^. This fluctuation 
is stretched by the cosmic expansion, and once the wavelength of fluctuation exits the 
horizon, it is frozen and treated classically. Quantum fluctuations continuously exit the 
horizon, and this mechanism creates the same structure of perturbations at every length 
scale (see e.g. [3]). 

The observations of WMAP [1] flnd a nearly scale invariant spectrum of primordial 
temperature fluctuations. It is often stated that WMAP conflrmed inflation, and the 
results expected from PLANCK satellite will narrow down possible models of inflation. In 
making such a statement, it seems that a particular mechanism [5] for generating CMB 
fluctuations is assumed, which is based essentially on slow-roll inflation [6]. Regarding the 
fact that theories of inflation have not been derived from fundamental theory of quantum 
gravity yet, we believe it is important to examine whether there are any issues that are 
overlooked in making predictions from inflation. 

In this paper, we point out that the temperature fluctuations of CMB can be generated 
by purely quantum effects, which is different from the mechanism usually assumed in the 
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slow-roll scenario. Our mechanism is based on the vacuum fluctuations of a large number 
of fields that are classically at rest at the bottom of their potential. The effect from each 
field is small, but we show that a sufficiently large number of fields from Kaluza-Klein (KK) 
modes or string excitations can produce an observable level of temperature fluctuations. 
The effect of these fields on the tensor perturbation is small. 

By comparing the temperature fluctuations obtained from our mechanism to the ob- 
served amplitude, we find that a theory with relatively low fundamental scale (i.e. string 
scale being 5 orders of magnitude lower than the 4D Planck scale) is favored. 

1.1 Temperature fluctuations 

To clarify the difference of our mechanism for generating temperature fluctuations form 
the one in slow-roll scenario, let us briefly review the latter [5]. 

In most of the currently studied models of inflation, it is assumed that a scalar field 
(inflaton) goes through a classical motion. In slow- roll inflation [6], inflaton rolls down the 
potential which is flat enough for vacuum energy dominates over the kinetic energy. In 
chaotic inflation [7], potential is generic but the friction due to the expansion makes the 
motion effectively slow. In other models such as N-flation [8], slowly moving classical field 
is effectively involved in certain sense. 

The equal value surface of the inflaton field provides a natural time slicing. Thus, 
fluctuations of inflaton 5ip can be reinterpreted as the fluctuation of time duration, or how 
much the universe has expanded: The slice of = is obtained by gauge transforma- 
tion, 5t = —Sip/ipci, and on that slice, fluctuation of the spatial curvature TZ = —H6t is 
generatecjl] (which can be translated into the gravitational potential $). In the slow-roll 
scenario, curvature perturbation is enhanced due to the slowness of the classical motion 
l/0ci- Using the fact that 6ip ~ H, and the slow-roll approximation, 0ci = V'/{3H), 
curvature perturbation is written as 7^ ~ {l/e){H/mp), where e = {VmpY / {SttVY is a 
slow-roll parameter characterizing the flatness of the inflaton potential (see e.g. 

During inflation, tensor perturbation is also generated [12]. In the linearized approxi- 
mation around an isotropic background, transverse-traceless tensor is decoupled from other 
fields. It satisfies massless equation of motion, and the amplitude is of order H/rrip, as is 
clear from dimensional analysis, where rrip is the Planck scale. Tensor perturbation will 

^Although this heuristic derivation gives the correct answer in the slow-roU limit, for a consistent 
analysis, one should use the gauge invariant variable defined in [9] which corresponds to the curvature 
perturbation on comoving hypersurfaces. This enables one to study the general cases. See e.g. [TU] . 
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produce B-mode polarization in the CMB. This is not observed at present, and the tensor 
to scalar ratio is bounded above by tt/s = 27^/7^^ < 0.2. This leads to an important 
conclusion that H is at least 5 order of magnitude smaller than mpi. 

We propose a mechanism for generating temperature fluctuations by vacuum fluctu- 
ations of the fields which are classically at rest. (Our mechanism is different from N- 
flation [S] in this sense.) Energy momentum tensor is quadratic in these fields, and their 
effect on the gravitational potential $ is neglected in the usual first order perturbation 
theory. Each field gives a small contribution of order {H/nipiY to $, but when there are 
many fields (typically of order ~ 10^^), this can sum up to an observable level. 

The fields with small mass compared to H do not oscillate during inflation, since 
the friction due to cosmic expansion overdamps the oscillation. These fields contribute to 
temperature fluctuations. When there are extra dimensions whose size L is large L ^ H~^, 
we have a large number of KK modes which contribute. The effect of these fields on tensor 
fluctuations is shown to be small. In our approach, the enhancement of scalar perturbation 
to tensor perturbation is due to the large number of fields that contribute to the former. 

In this paper we first compute fluctuations assuming the background is pure de Sitter, 
and later discuss the changes needed when Hubble is time dependent. Since the fluctuations 
originate from massive fields, the spectrum is tiled towards the UV (spectral index rig > 1), 
if Hubble parameter were constant. However, the spectral index is strongly dependent on 
the time-dependence of H. It can be lowered if Hubble decreases with time. We cannot 
know the dynamics of Hubble unless we know the origin of vacuum energy during inflation. 
In this paper, we do not make definitive statement, but we mention the possibility that 
quantum fluctuations of these fields (renormalized expectation value of energy momentum 
tensor) is the source of vacuum energy. 

Related work has been done by Nambu and Sasaki [13]. They computed correlation 
functions of energy-momentum tensor at the quadratic order in fluctuations, and related 
them to curvature perturbations. Their analysis is very similar to ours, but the setup 
and the interpretation are different. They consider a scalar field in an unstable potential 
< (with a suitable regularization). Their goal is to rederive density fluctuations 
in slow-roll inflation from purely quantum analysis without directly using the classical 
solution which rolls down the potential. On the other hand, we are considering fields in 
the stable potential rn? > 0, and studying their vacuum fluctuations. 
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1.2 Organization of this paper 

We will include descriptions of some known facts to make this paper self-contained and to 
clarify our assumptions. 

In Section 2, we review quantization in de Sitter background. In Section 3, we study 
Einstein equations and express gravitational potential $ in terms of matter fields. In 
Section 4, we obtain two-point functions of $. In Section 5, we find the CMB temperature 
fluctuations, and compare our formula with the observed amplitude to find typical value of 
parameters of fundamental theory. In Section 6, we study the spectral index, and discuss 
the effect of time-dependent Hubble constant. In Section 7, we consider non-Gaussianities. 
We compute three-point functions at the lowest order in the interaction, and we estimate 
the importance of interactions. In Section 8, we give a summary. In Appendix A, we 
perform the analysis of fluctuations including an infiaton field as an effective model for 
time-dependent Hubble. 

Part of the results of this paper has been reported in our previous publication |14j . 

2 Quantization in de Sitter space 

In this and the following two sections, we derive the formulas assuming the background is 
pure de Sitter space. We will discuss later what kind of changes are needed when Hubble is 
time dependent. We start by reviewing the calculation of correlation functions in de Sitter 
space, paying attention to the behavior in the small mass limit, which will be important 
for later applications. 

We will consider free fields, since we are mainly interested in weakly coupled theories. 
The magnitude of temperature fluctuations described in this paper depends on the number 
of fields which have masses smaller than the Hubble scale, but not on the details of the 
theory, so our conclusions will be valid even in the presence of interactions. We will discuss 
the effect of interaction in Sec. 7. 

The metric of de Sitter space is 

ds'^ = dt^ - a^{t)dx^, a{t) = H'^e^^ (-oo < t < oo) (2.1) 
= a\T) {dr' - dx^) , a(r) = -^-^-y, (-oo < r < 0) (2.2) 

where the conformal time r is defined by r = J dt/a{t) = — e~^*. 
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2.1 Scalar s 



Let us consider a free massive minimally-coupled scalar field, 



(2.3) 



It is convenient to define a rescaled field %(''", a^) = o(t)0(t, x) which has the standard ki- 
netic term. The equation of motion for the Fourier mode Xki'^) where x(t, ^) = J (Cjs Xfc(''")e*'^^ 
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2) -2 Xdr) = 0. 



(2.4) 



Canonical quantization condition is [x{t, x), x'(t, x')] = i6^{x — x'). We define the creation 
and annihilation operators at, by 



X[r, X) 



(27r)3/2 



2\k\ 



(2.5) 



where ^^(t) is the solution of (12.41) which is normalized as w^wi — u*4if: = 2i\k\. We take 
the solution which approaches u^{t) — > e^*'^'"^ at early time r — — oo, so that the choice 
of the vacuum reduces to the one for fiat spacetime in the short-distance limit. We will 
take Bunch-Davies vacuum, which is annihilated by a^'s in (12. 5p . throughout this paper. 
The explicit form of Uj^ir) is 



■\k\THl'^ 
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with 



(2.6) 



(2.7) 



Asymptotic behavior at the late times (in the super-horizon \k\/a ^ H limit) is given by 
^fe ~ ("1^1''")"'^^^ 5 i'^ terms of the original field, (p ~ {—t)^^'^~'^, as we can easily from 
the equation of motion (12. 4^ : The |/cp term drops out from the equation at late times, and 
the scaling w.r.t. time is independent of de Sitter symmetry tells us that the spatial 
(|A;|) dependence enters as a multiplicative factor with the same scaling dimension as the 
one for r (see e.g. [T5]). 
Fields with small mass, 
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do not oscillate in time. The friction due to the cosmic expansion overdamps the oscillation 
due to energy of massive field. 

We are interested in correlation functions, which are the expectation values taken with 
Bunch-Davies vacuum as in and out state. Two-point function at equal time is given by 

tf(r,f)^(r,?)) = -i-^ / (2.g) 

Substituting the late time expression for ^^(r), we get 

(0(r, x)0(r, x')) = H^C{rf){Ha\x - x'\)-^ (2.10) 

where 

7 = 3-2z/, (2.11) 
In the limit of small mass mH~^ ^ 1, we have 7 ~ ^m^H"^, and 

The coefficient C{'~f) diverges in the mH^^ — limit, but physical quantities such as the 
gravitational potential $ stays finite in this limit as we will see below. In the massless 
limit, the exponent 7 approaches zero, and the decay is slowest. We will see that fields 
with small mass (more precisely mH^^ < 10~^) mostly contribute to $. 
The energy-momentum tensor for a minimally coupled scalar is given by 

ST^u = S^d^(pd,(P - ^g^,{d'(j)d,(P - mV2)| . (2.14) 

This serves as the source for the gravitational fields. Let us look at the time-dependence 
of the 5Too component. Since (p ~ (— r)'''/^, the leading term of 5Too scales as 

5Too ~ i-ry-'. (2.15) 

We will see in the next section that this produces the gravitational potential $ ~ {—ry, 
which decays slowly when 7 ~ '^rn?H~'^ <^ 1. The fields that give important contributions 
are those which give 6Tqq ~ (^—-j-'^-'^+oim? h ^^^^ times. We can safely neglect the 

fields for which 5Tqq decay faster than this, when we compute $ in the late time limit. 
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So far we have considered minimally-coupled scalar. If there is coupling to the curva- 
ture, the action becomes 

S = \J d^x^/^ {d^<pd^^(t> - (^«' + iRW] ■ (2.16) 

where R is the scalar curvature of the background; ^(> 0) is a constant, which takes 
the value ^ = 1/6 for the conformally invariant coupling. For de Sitter space, we have 
R = 12H^. The curvature coupling effectively increases the mass by 

m^^m^ + 12H^^, (2.17) 

and changes z/ to 

u = - 12^ -m^H-^. (2.18) 
The late time behavior of such a field is 



0^(-r)l"VF^2r;^. (2.19) 

For the conformal scalar (^ = 1/6, m = 0), decays as ~ (— t)^- To have the exponent 
close to zero so that the field contributes to $, we need ^ = (minimal coupling) or close 
to zero, and mH~^ ^ 1. 

2.2 Vectors 

Massive vector field (Proca field) is described by the action 

S = I d'xV^g {^^g^^' g'^'"' F,,F,,,, - '^g>^^' A,A,}j . (2.20) 

Vector field arising from the KK reduction of a gauge field in higher dimension is such an 
example. 

The equation of motion in de Sitter space (in the conformal coordinates) is 

r]>'>''d^{d^,A, - d,A^,) + m\^A, = 0. (2.21) 
We act diy on this equation, and find a constraint, 

doAo + 2nA, = d,A, (2.22) 

where 

a' 1 

n = - = —. (2.23) 
a T 



To solve the equation of motion, we decompose Ai into the transverse and the longi- 
tudinal part, 

Ai = aP + dia (2.24) 

(T) 

where diA] =0. The transverse part satisfies 

{d^ - A + m2a2)Af ) = 0. (2.25) 
The component Aq satisfies the same equation after a rescaling by the scale factor, 

((9^- A + mV)(aAo) = 0. (2.26) 
The scalar function a is determined by 

A a = ^0^0 - 2—. (2.27) 
r 

These equations f l2.25p . (12.261) are equivalent to the equation of motion satisfied by 
^V^conf where v^conf is a scalar with conformal coupling, ni?H'^ — )■ m^H^ + 2. At late time 
and in the limit of small mass, the fields scale as 

Ao ~ ^conf ~ i-ry+oim^^-^) (2.28) 
aP ~ a^.o„f ~ (2.29) 

Also, from ([ISZD, we find a ~ (-r)°("^'^"'). 

The energy-momentum tensor for massive vector field is 

6T^, = F^pF/ - m'A.A, - ^g,, (^^p-^"" " m^A.A"^ . (2.30) 

To see the scaling of STqq at late times, let us look at its mass- dependent part. 



if - -^«VLnf - (2.32) 



2 

'^^oo"'^ = -^(^0^0 + AA,). (2.31) 
The leading time-dependence is given by 

2 

The energy momentum tensor scales in the way as if we had the scalar field </3conf • Whether 
the field contribute to $ or not depends on the effective mass in the equation of motion. 
Vector fields in (3+1) dimensions have conformal coupling (due to the conformal invariance 
in the massless limit), and decay faster than the minimally coupled scalars at late times, 
so they do not contribute to $. 



8 



2.3 Spinors 

The action of Dirac spinor is 

S = j d^x^^ ('^e^7'^^M -m)^>]. (2.33) 

In the background (12.21) . this is written as 

S = j d^x{{a"^^!) {zYdf, - ma) (0^/2^)} , (2.34) 

reflecting the fact that spinors are conformally invariant in the massless case. 
Dirac equation is 

{i-ff'd^ - ma)(a3/2^) = 0, (2.35) 

and the independent components (a^^^'^ = (i/j^^ip-) in certain representation of gamma 
matrices) satisfy 

dl-/\ + 5 U± = 0. (2.36) 



r 



Note that in the massless limit, fl2.36p is the equation satisfied by the conformal scalar 
a^pconi- Thus, the original field scales as \1/ ~ a'^^'^'^coni ~ (— 7-)3/2+o(m2H 2) ^^^^ times. 
Energy momentum tensor for spinors is (see e.g. [21]) 

5T^, = {^K^^D,)^ - (D(^^)7,)Vl/} . (2.37) 
The 5Tqq component scales as 

since 70 = e^o7'^ has one factor of a ~ (— r)~^, and decreases the power of r by 1. 
This 5Too is smaller than that for the massless minimally coupled scalar, so spinors do not 
contribute to $ at late times. 

2.4 The fields that are important at late times 

We have seen that minimally coupled scalar with mass rriH^^ ^ 1 decays most slowly, 

the late time limit. Coupling to the curvature effectively increases 
the mass, and fields such as conformal scalars decay faster. The fields whose independent 
components scale in the same way as minimally coupled scalar can contribute to $ in the 
late time limit. 

We can have small mass for the KK modes when extra dimensions are large enough 
L ^ H~^. Let us list possible origins of the fields which have minimal coupling. 
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• Massless minimally coupled scalars in higher dimensions. 

• The scalar fields from the KK reduction of gauge fields with indices along the internal 
directions: As long as the size of the extra dimension is stabilized independently 
of the scale factor for the 4D spacetime, these field do not have coupling to the 
curvature. Higher dimensional graviton with indices in the internal directions is also 
such an example. 

• Massive tensors (in 4D) from the KK reduction of higher dimensional gravitons: In 
the massless limit, the transverse mode satisfy the equation of motion equivalent to 
massless minimally coupled scalar (see (13.251) below), thus contributes at late times. 

Whether the first type of fields exist or not may depend on the theory, but the second 
and the third (one-form gauge fields and gravitons) will exist in fundamental theories 
in general. In the following, we will not ask how many of these fields exist. We will 
ignore order 1 factor coming from this multiplicity, since this is much smaller than the 
huge multiplicity of the KK modes for each field. In the explicit analysis, we will take 
minimally coupled scalar fields. Other fields can be studied in the similar manner by 
considering the independent components which satisfy scalar-type equations of motion, as 
long as we are considering vacuum fiuctuations of these fields. 



3 Einstein equations 

We now study Einstein equations. Einstein equations are constraint equations which allow 
us to write the gauge invariant metric fiuctuations (such as gravitational potentials $ and 
\I^) in terms of matter fields. The metric fiuctuations are decomposed into scalar, vector, 
and tensor modes, each of which can be studied separately. Tensor mode is the part 
which is transverse-traceless in the spatial directions, vector modes are those which are 
divergenceless, and scalar modes are those which can be written as derivatives of scalar 
functions. We follow the notation of [1]. 
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3.1 Scalar fluctuations 

The scalar part of the (0,0), (0, i), (z, j) components of Einstein equations are given, 
respectively, by 

- 3H{^' + H^) = AttGSToo, (3.1) 
(*' + ?^$),i = 47rG5Tof , (3.2) 



+ H{2^ + $)' + {2H' + ?/^)$ + y ($ - *) 



5..-^(<^'-*),..=4vrG5rf . 



(3.3) 



We take the background spacetime to be pure de Sitter space (a = —H^^/r, Ti = — 1/r). 
The l.h.s. are the Einstein tensor expanded to the 1st order in metric fluctuations. $ and 
are the two gauge invariant variables constructed from the scalar components. In the 
longitudinal gauge, they are given by 

dslg. = a^{{l + 2^)dT^ - (1 - 2-(I/)Sijdx'dx^} . (3.4) 

On the r.h.s., we take the energy momentum tensor which is quadratic in the matter 
fields. We consider minimally-coupled scalars here. We assume there are many free scalar 
fields. The energy momentum tensor is a sum over their contributions. 



(3.5) 



For brevity, the label on the field is suppressed, and the sum is understood to be over the 
species. The fields have vanishing classical background, and are gauge invariant. Each 
component of 6T^i, is given by 

SToo = + 9,09,0 + mV02}, (3.6) 

STo, = ^{0'9,0}, (3.7) 



STi, = J2 {^^4>^J4> + ^6ij (0'2 - 9,09,0 - m2a202) | _ 



(3.^ 



The superscript (S) in ( 13. 2p . ( 13. 3p denotes the scalar part. Recall that STqh and 6Tij can 
be decomposed as 

SToi = diS + Ui, (3.9) 
5Tij = didjS - -6ij A s + diVj + djVi + tij + fSij, (3.10) 

11 



where Ui and Vi are transverse vectors diVi = diUi = 0, and tij is a transverse traceless 
(TT) tensor, ditij = tu = 0. By the scalar part, we mean the part involving s in (13.91) . and 
the part involving s and / in (I3.10p . 

We can find s and s by taking divergence and applying inverse Laplacian, 

s = ^^dkdi{STki - ^SkiSTmm) = ■^^didj{di(t)dj(l) - ^dk(l)dk(t)). (3.12) 

Using Einstein equations (I3.1l) - (l3.3p . we can solve for $ and \i/ in terms of 0. First, 
from the traceless part of (13. 3p . we find 

$ - ^ = -SttCs. (3.13) 

Using this in (13. 2p . 

^' + H^ = 8nG |-/ + ^9,(0'9,0)| . (3.14) 

The last term is the part that we would get if we had $ = 

Let us solve (I3.14p by substituting the late time asymptotics of (p on the r.h.s. This 
is a valid procedure, since we are interested in the correlation functions of $ in the late 
time limit, and $ only appears as external lines. Special care is needed if the leading 
term (which has the lowest scaling dimension) is degenerate with another term, which can 
happen at certain values of the parameter; we will comment on this point when necessary. 

In the late time limit, the time- and space- dependence of the field factorizes, 

0(r,x) = (-r)^0(x), (3.15) 

where 7 is defined in (12. lip . Time-dependence of $ is found from the time-dependence of 
the r.h.s. of ( Km . 

r(l$)' ^ (-r)^~i $ ~ {-ry. (3.16) 

r 

Thus, $ at late times can be written as 

$(r,x) = (-r)^$(x). (3.17) 
Time derivative is given by $' = \E'' = ^\E', 0' = ^(p. From (I3.14p . we get 

$ = 471^^ |_^a,a,(a,0a,-0) + ^d,^d,^ + ^A , (3.18) 



12 



where we have used di{(j)di(j)) = A0^/2 to rewrite the last term. This solution is consistent 
with all the other components of Einstein equations. 

The expression (13.181) diverges at 7 = 1. At this special value, fl3.14p cannot be solved 
with the naive ansatz $ ~ (—r)'^, since the l.h.s. vanishes. In this case, we can solve the 
equation by setting $ ~ (— t) log(— r). 

Also note that there is always a freedom of adding a term which has time-dependence 
~ (— t)^ to the solution of (13.181) . but we can eliminate this piece by requiring that the 
solution does not blows up in the early time limit. 

3.2 Vector fluctuations 

The vector modes are the following part of the metric fluctuations, 

ds"^ = [dr^ + 2Sidx'dT - - Fij - Fj,,)dx'dx^] , (3.19) 

where diSi = diFi = 0, and ,i denotes derivative w.r.t. x\ There is a gauge invariant 
combination, 

Vi = Si- F[. (3.20) 
The vector part of the (0, i) and (z, j) components of the Einstein equations are 

A 1/^ = 167rG5To^f \ (3.21) 
{Vid + Vj.i)' + 2^(^*J + V,-,.) = IQtiG5T'>J\ (3.22) 

where the superscript [V) denotes the vector part, which are the part of dT^i and 5Tij 
which involves Ui and Wj, as defined in (13. 9p and (13.101) . Ui is given by subtracting the 
scalar part from 5Toi, 

dT^'^^ = Ui = cp'd4 - ^dM<P'dk4>)- (3.23) 

The leading term of Ui at late times is smaller than it naively looks. Recall that 
~ (_7-)7/20(i _^ 0{t^)). We can see that the order {-t)^~^ term of the r.h.s. of (1323]) 
vanishes by using 0' = ^0 and (pdicj) = di{(fP')/2. Thus the leading term of the r.h.s. of 
(I3.2ip scales as (— r)'^"*"^, which implies Vi ~ (— r)'^+^. This behavior is consistent with the 
equation (13.221) . 

Since Vi decays at least as (— r)^, we conclude that the vector perturbation produced 
by the matter fields (j) is negligible at late times. 
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3.3 Tensor fluctuations 



The transverse-traceless (TT) tensor fluctuation hij {y^hij=h\ = 0) is defined by 

ds^ = [dr^ - {5ij - hij)dx'dx^] . (3.24) 
It is sourced by the TT part of energy-momentum tensor, 

h'l^ + 2m\^ - Ahij = 87rG6Tp. (3.25) 

The general solution to this equation is given by the solution hf^^ for the homogeneous 
equation on top of a particular solution hlj which depends on dT^ ' . 

The homogeneous equation is equivalent to massless scalar equation of motion. Its 
solution hf^ is the usual gravitational wave, which scales logarithmically in space and 
time. This has the scale invariant spectrum with the amplitude H/nipi. 

The time-dependence of h^^^ is determined by (13.251) to be h^^^ ~ (— r)'^"'"^, since ST^p ~ 
(— r)'^. Care is needed when 7 = 1. In this case, (— t)'^"''^ = {—t)^ is degenerate with the 
(decaying) solution of the homogeneous equation, and fl3.25p cannot be solved with this 
ansatz. In this case we have to take h\^^ ~ (— r)'^ log(— r). In any case, /i-j^ decays at late 

(T) 

times, and the effect of 6T^- for the tensor fluctuations is negligible at late times. 



4 Correlation functions 

Having expressed $ in terms of (p, it is straightforward to compute correlation functions 
of Let us compute the two-point function. 
We decompose $ in (I3.18P into two pieces, 

$ = $o + $i (4.1) 
1-7 



$1 = AnG^{—d,d,{di<pdj^)--id,^d,^)], (4.3) 



where $0 is the part which we would get when $ = and $1 is the part which depends 

on s defined in f l3.12p . 

The ($o'^o) correlator is just a product of two propagators, 

2 

(0(r,x)0(r,xO)' 



($o(r,x)$o(r,x')) = '^inGYj^^Y^) 

= 2(vrGi72)2 5^(^^yC2(7)(-r)2^|a:-a:'|-'^ (4.4) 
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where as in the last section, the sum is taken over all species of 0. When the mass of the 
fields are small (7 ^ 1), we get 

{Mr,x)Mr, x')) ~ ^(G/j2)^5^(-r)2^|x-a;'|-'". (4.5) 

The contribution from the fields with 7 ^ 1 is finite, since the two factors of I/7 from the 
propagator are canceled by the two factors of 7 from fl4.2p . 

The other parts of the correlator can be computed by using the formulas such as 

(S'i = i;) 

SAr^ - A|A.-(A.2) '--r-":;--n (4.6) 



3 I rf rf^ I 11 />"» rpf I A~)~2 I rr* rff I A~[~4 

1 A(A-l) 
which are valid up to possible contact terms. The cross term ($i$o) is 

(<l>i(r,x)$o(r,x')) = -K^GH^Y^ir^^ [|^9,a,(a,0(r, x)0(r, x'))(9,0(r, x)0(r, x')) 

--^{di<p{T, x)(f){T, x')){di(j){T, x)(j){T, X')) 

= -4(vrGi7^f j:(^y(^)c^(7)(-r)^1x-x'|-^. (4.8) 

In the following, we will find that the fields that mostly contribute to the CMB tem- 
perature fluctuations are the ones with mH^^ < 10"^, so let us study the 7 <^ 1 behavior 
here. In this limit, (14.80 is smaller than (I4.5P by a factor of 7. 

The part 

2 ^2^272 + 47-3) 



($,(r,x)$,(r,a:0) = 16(vrGif^)^ E (l^:^) 



L)(7 + 3)(27-l)(27+l) 
xC2(7)(-r)2^|x-x'|-'^. (4.9) 

is smaller than (14.51) by a factor of 7^. Thus, we can consider only the ($0^0) psj:t when 
7 < 1. 

4.1 Summing up KK modes 

Let us perform the summation over the massive fields </>, assuming they are Kaluza-Klein 
modes from the compactification of extra dimensions. For definiteness, let us assume there 
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are D dimensions which are compactified on a torus with the periodicity L (assumed 
to be the same for all directions for simplicity) being large compared to the inverse Hubble 
of inflation, L ^ H~^. We assume the internal directions other than these D dimensions 
are compactified on a space with the string scale size. 

The mass of a KK mode with {ng} units of momentum on is 

a=l 

When the level is sufficiently dense, the density of states in the mass interval dm around 
m is given by 

SD^i\n\'^~^d\n\ = 5z)_i(L/27r)^m^"Mm, (4.11) 

where Sd-i = 2n^^'^ /T(D /2) is the volume of the D — 1 dimensional unit sphere. This 
relation states that the number of states is proportional to the phase space volume (the 
volume element of the KK momentum space times the volume of the internal space). 
Using (14.111) . we convert the sum in the ($$) correlator to an integral, 

/ dmm^-\Ha\x-x'\)-^'', (4.12) 

mpij Jo 

where cd = S'£)_i/(4(27r)^"'"^). We have used the expression (14.51) for the correlator in the 
7 <^ 1 limit. The upper limit of the integration should be ~ as long as we are 
working in Einstein gravity, so that the field </> are the ones which do not oscillate. 

However, if string scale is less than Hubble scale, < H, string states also have to 
be taken into account. In this case, we expect that the sum over the mass is effectively 
cut off at rric ~ nig for the following reason. Let us assume the two-point function of $ 
comes from the one-loop diagram in string theory (Fig[T]). String theory can be regarded 
as a field theory with infinitely many fields, except that one-loop amplitude effectively 
has UV cutoff due to modular invariance. The integral over the moduli r is restricted 
to the fundamental domain (Fig. [2l), and the Schwinger proper time (Im(r)) is cut off 
at string scale. There is no physical meaning to time interval shorter than string scale, 
or oscillations much higher than string scale. This means the internal states in the loop 
which has mass much larger than string scale do not have physical effect. This argument 
is based on the perturbative string theory in fiat spacetime, and it is not clear whether it 
is valid in an arbitrarily curved background, but we believe this is a reasonable estimat^. 

^In fact, the precise value of the upper hmit of integration is not very important in the parameter 
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One-loop in field theory String torus diagram 



Figure 1: One-loop diagram: One- loop diagram in field theory corresponds to the torus 
diagram in string theory. String theory can be regarded as a field theory with an infinite 
number of fields, which are the Fourier modes on the world sheet spatial direction. 

5 CMB temperature fluctuations 
5.1 Generation of adiabatic fluctuation 

The formula fl4.12p . which was obtained with constant if, is approximately valid during 
inflation. H decreases towards the end of inflation. When H becomes less than the mass 
of a field 0, the field will undergo classical oscillation, and decay into radiations and 
stable particles. The energy density from this process produces curvature perturbation. 
This is similar to what happens in the "curvaton scenario" [16], in which the curvaton 
field (different from infiaton, and usually assumed to be one field), produces curvature 
perturbation. The situation that we are considering has similarities and differences with 
the situation usually considered in curvaton scenario. Even though it is helpful to have 
curvaton mechanism in mind, we emphasize the following specific assumptions that we 
make. 

We assume that the KK modes and string states decay sufficiently fast so that they do 
not interfere with the standard big bang nucleosynthesis. The decay products are assumed 
to be in local thermal equilibrium, and can be treated as a single fluid which dominates 
the energy density of the universe. We do not assume low-energy (such as TeV scale) 
supersymmetry. In such a case, it will be generically expected that the decay products 
interact among themselves, reaching thermal equilibrium. This is the reheating mechanism 

region of interest. The conclusion that nisH^^ ^ 0.1 is favored does not change even if we take the upper 
Uniit to be {3/2)H instead of nis- 
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Figure 2: Fundamental region: Integration over the moduli of the torus should be restricted 
to the fundamental region, since other regions are physically equivalent to this region due 
to modular invariance. Im(r), which corresponds to the Schwinger parameter, is cut off 
at the distance of order string scale. 

in our model. The particles that have conserved quantum numbers at present, such as 
baryons, cold dark matter, etc. will be produced after the above thermal equilibrium has 
been established. In this case, fluctuations of the density of these species are determined 
solely by the local temperature, and obey the adiabaticity relation (see e.g. discussion 
in [19]). In the context of curvaton, this corresponds to the case where matter (such as 
cold dark matter and baryons) are produced after curvaton decays, in which case there 
is no isocurvature fluctuations [18]. Even though the mass scale of KK modes is lower 
than the Hubble of inflation, it will be much higher than the scale of Standard Model of 
particle physics. In such a case, we do not have reason to expect that the fluctuations 
of particular species of KK modes to be directly related to those of cold dark matter or 
baryons. Thus, we do not expect there to be significant isocurvature perturbation (which 
is strongly constrained by observation; see e.g. Section 3.6 of [IT]) in our model. 

The amplitude of the energy density fluctuations produced by the above process is 
determined by the quantum fluctuations of 0's during inflation. We assume the transition 
from quantum fluctuations to classical oscillation and the decay of these fields occurs 
quickly and more or less simultaneously for all the fields. (This is a simplification to 
make the argument simple; we leave more general analysis to future work.) Under this 
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assumption, we identify $ evaluated at the end of inflation as the one sourced by the 
thermal fluid described above. Once the universe is in thermal equilibrium with a single 
fluid dominating the density, the superhorizon mode of $ basically remains constant (see 
e.g. [1]). It changes only by order 1 factor at the transition between matter and radiation 
domination, but ignoring this factor, $ at the end of inflation is directly related to its 
value at the recombination. 

One may worry that this mechanism produces anisotropy, since there are fields with 
non-zero spin (such as KK modes of gravitons) whose components separately undergo 
classical oscillation. However, note that the inhomogeneities 6p/p that results from the 
fluctuations of many fields scales as 1/ -y/iV as the number of fields increases. This is 
because the possible classical density p will be proportional to A^, but the fluctuations 5p 
(or more precisely, the root of the square expectation value a/ will be of order ^/N 
if the fields fluctuate independently. Therefore, total anisotropy from the fluctuations of 
many fields scales as 1/ \/N and is kept small. This type of suppression of anisotropy due 
to a large number of fields appears also in the context of "vector inflation" |20] . 

5.2 Amplitude of the CMB temperature fluctuations 

Having stated our assumptions which lead to the identification of $ at the end of inflation 
with $ at recombination, let us now study the CMB temperature fluctuations 5T /T. 

Temperature fluctuation of CMB is related to $ at recombination (at redshift z ~ 1100) 
by 5T /T = — <l>/3. The angle 6 on the sky corresponds to the distance dr = 2RrSm{6/2), 
where Rj. is the radius of the surface of last scattering. This is of order the inverse of the 
present Hubble parameter ~ Hq^. The modes outside the horizon at recombination 
correspond to the angle 3° < ^ (or angular momentum I < 60). These modes have been 
outside the horizon since the inflation, and $ is frozen (remain constant at fixed comoving 
distance). The distance dr corresponds to the distance 

{ae/ar)dr = 2Rsm{9/2) = a^lx — x'\ (5.1) 

at the end of inflation. The radius R = {ae/aj.)Rr will depend on the scale of inflation. We 
will take the standard estimate RH ~ 10^^ ~ e^'' in the following, which amounts to the 
assumption that the reheating temperature is not much lower than the grand unification 
scale. 
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The angular power spectrum Ci is defined by 



:^W?^(0)) = £(2/ + l)QP,{cose). (5.2) 



T 

1=1 



We will focus on the superhorizon modes. To find amplitude of these modes, we expand 
the coordinates in the correlation function fl4.12p (around 6' = tt) as 

{Ha\x-x'\y^^ ~ (2i?iy)-2^(l -27log(sin(^/2)), (5.3) 

and recall that — 2 log(sin(^/2)) is !/(/(/ + 1)) in harmonic space. From fl4.12p . we find the 
square amplitude 6^ = l{l + l)Ci, 



(5.4) 



Note that the integrand is strongly suppressed at m > due to the factor (2i?if )~t™^^ ^ . 

In (15.41) . we have taken the upper limit to be ttlc = m^. This is a good approximation even 
when rris > H, since the integral has little contribution from the region near the upper 
limit. We have defined 



Md{0 = dte-^'\^+\ ^0 = \og{2RH). 



(5.5) 



D+2 

To see the qualitative behavior of it would be helpful to note A^d(Co) ~ Co ^ 
when Co > 1- In this limit, we have 6"^ ~ {—)^{LH)^ {\og{2RH))~^ up to constant 
factors. 6^ is enhanced when extra dimensions are large, (LH)^ ^ 1, since many fields 
contribute to it. 6'^ becomes small if \og{2RH) were larger due to the decrease of massive 
wave function at large separation. 

5.3 Comparison with the data 

We will now use observational data [4], 

5t ~ 2.6 X 10-^ n/, < 0.22, (5.6) 



to constrain the parameters in our model. This implies ^ = y ^^x'^t/s ^ 0.81 x 10 ^. 
Let us first assume this inequality is saturated. Then the amplitude (15.41) provides the 
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relation between the two parameters and L, or equivalently, between and the string 
couphng gs, since L is written as (Lm^,)^ = Sn^gglmpi/ml). 

Figs. [3] and m show (Lnipi) and gs as functions of nis/H, respectively. It is easier to 
have weak coupling with small D, while it is easier to keep L not too large with large D. 
Typical values that are consistent with (15. 6 p would be: 

{D = 2, m jH = 0.2, Lm^i = 10^^^ g^ = 3}, (5.7) 
{D = 3, mjH = 0.2, LrUpi = 10^°, g^ = 5}, (5.8) 
{D = 4,mjH = 0.1, Lrripi = lO^g^ = 7}. (5.9) 

The number of the fields that participate in 6T/T is roughly ~ [Lnis)^ . For the above 
choice of parameters, 10^"^ ^ N < 10^^. 




0.2 0.4 0.6 0.8 1.0 

nis/H 



Figure 3: logj^o(-^"^pO a function of nis/H, with 6t = 2.6 x 10 ^, rt/s = 0.22, RH ~ 10' 



6 Time-dependent Hubble 

So far we have studied the fluctuations generated during inflation, assuming the back- 
ground is pure de Sitter. Let us now consider the effect of time-dependent Hubble. 
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Figure 4: gs as a function of nis/H, with 5t = 2.6 x 10 ^, rt/s = 0.22, ~ 10 
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6.1 Spectral index 

Since 6T/T originates from the fluctuations of massive fields, the spectrum is stronger in 
the UV. The spectral index Ug is shghtly larger than 1, 

d ^ /6T, ^^6T_ 



1 - 



d\og{Ha\x — x'\) 
4m2 MoiCo) 



'6T 



log(^(r,f)^(r, x') 



(6.1) 



3H^MD-2{Coy 

This is in the range 1 < n, < 1.02 when D = 2, and 1 < < 1.05 for D < 6 (See 
Figure [5]) . 

The above values are obtained by assuming the Hubble is constant. However, Ug is 
sensitive to the time-dependence of H. As long as the change of Hubble is adiabatic 
\H\/H'^ <^ 1, it would be reasonable to assume the amplitudes are determined in terms 
of H at the time of horizon crossing. We will have to replace the prefactor {H/mpiY in 
the ($$) two-point function (14.121) to [H {t\ior) / ^piY where thor is the time of the horizon 
exit for the scale of interest (e"-^*''™ ~ |x — a;'|). Since the long wave length mode exits 
the horizon early, the amplitude is lifted in the infrared. The spectral index is lowered by 
0.5 if there is time-dependence of order H / ~ —0.01. 

It would be necessary to understand the origin of vacuum energy during inflation to 
understand its time-dependence. In this paper, we cannot make a definitive statement. 
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Figure 5: Spectral index function of rris/H. 
but we would like to mention a possible origin of vacuum energy in the next subsection. 
6.2 Possible origin of vacuum energy 

In most inflationary models, the presence of vacuum energy (or a nearly fiat infiaton 
potential) is simply assumed and its origin is not clear. Also, in the recent constructions 
of de Sitter vacua in string theory, the mechanism for uplifting from supersymmetric 
vacua to de Sitter vacua is not fully understood. In the study of low energy effective 
action of string theory, it has been very difficult to find de Sitter vacua in a controllable 
approximation. (See e.g. [22] for a recent discussion.) 

It might be necessary to understand vacuum fluctuations of the fields to find de Sitter 
vacua. In string compactification models with large internal space [23], which is believed 
to be realizable generically, there are many light KK fields. Quantum fluctuations of these 
fields might be an important source of vacuum energy. 

With this motivation in mind, in this subsection, we discuss a possible dynamical 
scenario in which vacuum fluctuations and Hubble are determined in a self-consistent 
manner. 

Consider the expectation value (one-point function) of the energy-momentum tensor. 
This quantity is UV divergent, and we will renormalize it so that it vanishes in the fiat 
background. Because of de Sitter symmetry, the expectation value is proportional to g^j,^. 
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The renormalized expectation value of energy-momentum tensor of a scalar field in de 
Sitter background is given by (see (6.183) of [2T]) 

(T,.)ren = 1^ p jm^ + (e - |^(^ + + -y)- log(12m2i?-i)| 

_ ^2.^ _ \^ _ 1 2^ _ _ 1)2^2 ^^2l 2) 

6^ 18 2^^ 6^ 2160 J ^ ^ 

where V'(^) = r'(2;)/r(z), and v is the order of Hankel function as described in Section 2. 
(^Atj/)ren is rcuormalized by subtracting the divergent piece in the flat background. When 
there are scalar fields (assuming massless minimally coupled), we have 
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It would be possible that de Sitter space during inflation is a self-consistent solution of 



^i>.v - 'i^Qp.vR = -S7rG{T^^)ren- (6.4) 
The condition that both sides of the equation balances implies 

H'^N^, — (6.5) 

In fact, our scenario is not consistent with this equation as it is. The value of 
which produces the observed level of temperature fluctuations is at least ~ 10^^ with 
H/nipi ~ 10~^, and is too large by a factor of 10^ for (16.41) to be satisfled. However, 
this is not a contradiction. We do not expect that Einstein equation (16. 4p is applicable, 
since our preferred value of string coupling is nis/H ~ 0.1, and the left hand side will be 
corrected by string {a') effects, which are not negligible when rris/H < 1. 

Presumably, the self-consistent de Sitter solution of (16.41) is an unstable solution, and 
small fluctuation of H will drive the background to flat space, which is another solution of 
this equation. To study time-dependence, we will have to compute the expectation value 
(T^^)ren in the background with H ^ Q. This with (16.41) will tell us the evolution of Hubble. 
We will leave this analysis as an important open question. 

It is not clear whether the dynamics of quantum expectation value (T^,y)rcn is similar 
to the dynamics of inflaton, but let us assume it is for the moment. In Appendix A, we 
perform the analysis of fluctuations including inflaton fluctuations 5ip. We take 5ip to be 
of the same order as The gravitational potential $ has a term induced by 5ip in 

addition to the term from the matter flelds that we have studied. The relative importance 
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of inflaton and the matter fields depends on the details, such as the slope of the inflaton 
potential and the time between horizon crossing and the end of inflation. The effect of the 
matter fields will be important unless the slope is fine tuned to a small value. 



7 Non-gaussianities 

Non-gaussianities appear in a characteristic manner in our mechanism. We first describe 
the calculation ignoring interactions among the matter fields 0. We then remark that the 
magnitude of non-gaussianities is controlled by the coupling constant in higher dimension. 

The three-point function of $ is given by the triangle diagram where each pair of points 
is connected by (00) (see Figure 7), 

($(r,x)<l>(r,y)$(r,z)) (7.1) 

1 f 
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We define the non-linearity parameter /nl by a local replacement, $ — )■ $g + /nl'^'^ 
with a gaussian field Let us consider three-points at super horizon separation, and 
estimate /nl by expansing the coordinates as in (15. 3p . The local form of /nl is enough to 
characterize the magnitude of non-gaussianity in this approximation. From (17.11) . 

This is proportional to rt/s, and further suppressed by the other factors (See Fig [6]). 
For rt/s = 0.22, we have /nl < 10"^ 

The reason for the smallness of non-gaussianity is that ($$) is roughly proportional to 
the number of fields A^, and ($$$) is also proportional to in our setting. This makes 
the non-gaussianity small /nl ~ ($$$)/($$)^ ~ N^^ in the large limit. This is in 
contrast to the curvaton case [16], where non-gaussianity is necessarily large if curvaton is 
the only source of curvature fluctuations [T8] . 

When there are interactions among 0's, we will have higher loop diagrams, such as the 
ones in which a propagator traverses two sides of the triangle. (See Figure 7.) Even though 
there are many fields, interactions do not necessarily make /nl huge. Since the fields 
are KK modes or string states, there will be conserved quantities such as momentum in 
the internal space or the excitation number of strings. The third field in the right diagram 
in Figure 7 is determined by the first two, and this diagram will be of order X^N'^. The 
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Figure 6: Non-gaussianity /nl as a function of nis/H. 

natural magnitude of the coupling A will be \q/\/V where Aq is the coupling constant in 
higher dimension, and V is the volume of internal space. This is because the canonically 
normalized KK fields has a factor 1 / \/V relative to the higher dimensional field, and 
the interaction Xcf)^ has one more cf) than the kinetic term. Since ^ ~ as mentioned in 
Section 4.1, the factor A^A^^ associated to the right diagram of Figure 7 is just given by 
Aq. (This is equivalent to saying that the correlation functions are computed by Feynman 
rules in the higher dimensions with the diagrams of the type of Figure 7.) 

If Ao takes finite but small enough value that perturbation theory is applicable, non- 
gaussianities will be given by the diagrams such as the ones in Figure 7. The shape 
(momentum dependence) [25] of this type of non-gaussianity will be different from those 
arising from the usual slow-roll inflation [26]. The difference will be seen by studying the 
three-point correlations at subhorizon separations. We will leave this analysis to future 
study. 

8 Conclusions 

We have shown that the CMB temperature fluctuations (adiabatic perturbation) can be 
generated by the purely quantum effects of fields which are classically at rest. When there 
are a large number of fields, this can produce observable level of fluctuations. Tensor 
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Figure 7: Diagrams for the three-point function of $. Left paneh Ignoring the interaction 
among the fields 0, this is of order (the number of fields). Right panel: Higher loop 
diagram, where 0i and 02 are independent, but 03 will be determined from the conservation 
law (such as KK momentum conservation). This is of order X^N"^. 

fluctuations are hardly affect by this effect, and will remain of order H/rripi. In our 
mechanism, the enhancement of scalar fluctuations relative to tensor fluctuation is due to 
the large number of fields involved, and not due to the smallness of the slow-roll parameter 
as in the usual slow-roll scenario. 

When the size of the extra dimensions are large compared to the inverse Hubble during 
inflation, we have a large number of Kaluza-Klein modes which contribute to this effect. 
String excited modes also contribute if < H. We compare our results with observed 
amplitude, and find that rris/H ~ 0.1 is preferred. The size of extra dimensions is typically 
of order 10^ GeV~^ or smaller. 

There have been models of inflation based on TeV scale supersymmetry (see [28] for 
a review). Inflation and reheating have been studied in explicit string compactification 
in [27]: Supersymmetry is broken by the hidden-sector branes wrapped around internal 
cycles, and inflaton is given by closed string moduli. In that case, it is important to make 
sure that the decay of inflaton reheats the visible sector dominantly and not the hidden 
sectors significantly, which has been checked in [27| . This is necessary to avoid cosmological 
problems, such as the generation of large isocurvature perturbations. 

In our case, we do not assume low-energy supersymmetry. In this case, it will be 
generically the case that massive fields, such as KK modes, decay and reach thermal 
equilibrium. This will occur well before the standard big bang nucleosynthesis begins. In 
this situation, there will be no isocurvature perturbations, as shown in |19] . 
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We performed our analysis assuming the extra dimensions are compactified on a torus 
. We believe this captures the qualitative features of quantum effects of KK states in 
the general compactifications with L ^ H~^, whenever the multiplicity of the KK modes 
is similar to that for . There have been studies on string compactifications which 
realizes supersymmetric vacua with all moduli fixed. It is argued that the large-volume 
compactification is generically achievable in the construction of [23] • Understanding of the 
mechanism for uplifting to de Sitter vacua or realizing inflation is at a more qualitative 
level at present. Brane anti-brane pair [29] will be a candidate for such a mechanism. We 
expect our mechanism for generating CMB temperature fluctuations should be relevant in 
these contexts. 

Non-gaussianity in our mechanism is given by triangle diagrams, with possible cor- 
rections. The magnitude is controlled by the coupling constant in higher dimension. It 
would be possible in principle to distinguish our mechanism from others by the precise 
measurement of non-gaussianities. 

The main purpose of this paper is to study fluctuations without asking the origin of 
vacuum energy during inflation. But as we mentioned in Section 6, it would be possible 
that vacuum fluctuations (the renormalized expectation value of energy-momentum tensor) 
of a large number of fields are the source of vacuum energy. It is important to understand 
the dynamics of this vacuum energy [30]. We have included the analysis of fluctuations in 
the background with time- dependent Hubble by introducing inflaton in the appendix, but 
it is not clear to what extent this captures the dynamics of quantum vacuum energy. 
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Appendix 



A Perturbations including inflaton 

In this appendix, we assume the dynamics of time- dependent Hubble is effectively described 
by an inflaton field and study fluctuations including the inflaton fluctuation 5ip in 
addition to matter fields 0. 

We assume inflaton and matter are not directly coupled with each other. Thus, the 
equation of motion for is the same as before, except that Hubble is now time-dependent 
and % and a are not those for pure de Sitter, 

[dl + 2'Hdr - A + aW] = 0. (A.l) 

This does not couple to fluctuations of other fields, so the quantization of can be done 
at first. 

The action of inflaton is 

S = j d'x^g {^-g^^d.^d^v - V{v)^ . (A.2) 

We decompose ip into classical part (which is homogeneous in space) and fluctuations, 

(/9 = V5o(t) + 5v9(r,f). (A. 3) 

Classical part satisfies the equation of motion 

dlv^ + 2mv^o + = 0. (A.4) 
To find the equation of motion for dip, it is convenient to take the longitudinal gauge, 

ds^ = [(1 + 2$)dr2 - (1 - 2^)dx'\ , (A.5) 
and expand the equation of motion, 

{^gg'^'^d.^) + = 0, (A.6) 

to the first order in 5^9, $, \E'. Then, we get 

{dl + mdr - A + aV^^] bLp - (v9[,)($' + 3^0 + 2aV^$ = 0. (A.7) 
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Energy-momentum tensor for (/? is 



The classical part (zeroth order in 5ip) is 



2a2 
1 

2^ 



(K)' + 2V^) 



y(v,o)o^ 

rp{tf,fi)i ^ 

The linear part in 5ip is 



(A.8) 

(A.9) 
(A.IO) 

(A.ll) 
(A.12) 
(A.13) 



Note that 6T^'^^^j does not have off-diagonal components. 

To write Einstein equations, we include the linear terms in 6ip and the quadratic terms 
in in the energy-momentum tensor. The (0,0), (0, i) and {i, j) components of Einstein 
equations are as follows: 



(A.14) 
(A.15) 



^" + H(2^ + $)' + (27/' + H^)<l> + — ($ - ^) 



5^^. __($_^),.. 
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(A.16) 



where the summation is taken over the species of matter fields as in the main text. 

Let us study the leading behavior in the super- horizon limit following [1]. We assume 
inflaton is classically slow-rolling. 



m^o + v;^ = 0. 



(A.17) 
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We consider inflaton e.o.m. ( IA.7I) and the (0, i) component of Einstein equation (1A.15|) . 
In terms of physical time t (and in the slow- roll limit), 

3Hdip + V^^6^ + 2V^^ = 0, (A.18) 



(A.19) 



where we have ignored the terms in the energy momentum tensor of which are small in 
the small mass limit, 7 = |m^iJ~^ ^ 1 (we ignore the term s' in f l3.14p ). 
These equations imply 

|«)-Ef/W---0, (A,20) 
where we have defined /(0) by 

19,(09.0) = -2|:/(0)e-^^*. (A.21) 

That is, 

/(0)e-7^^* ~ ^0' (A.22) 

in the late time limit. From (]A.20p . 

5y, = ^(c-lj]/(0)e-^^*) (A.23) 

with a constant C. We fix C so that the amplitude of 5ip is H at the horizon exit (since 
5ip is essentially a massless scalar inside the horizon), 

^{//.(^) +l^/(0)(e-'r^*.-e-^O}, (A.24) 

where the star denotes the quantities evaluated at the horizon exit. 

In the usual slow-roll inflation, there is only the first term. The factor {^i/^ generally 
grows towards the end of inflation, and it is assumed to be of order nip at the end of 
inflation. The second term represents the effect of matter fields to the evolution of inflaton 
fluctuation. 



5ip 



In terms of $, 



'^^j:'-^me-^'' (A.25) 

+ (A-26) 



H ^ 2 
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The first term is the part induced from the inflaton fluctuation (1A.24P through the usual 
mechanism. The second term is the effect of matter (agrees with the formula that we have 
obtained), which exists even if there is no inflaton. 

This expression will be valid until the end of inflation. After inflation, $ will be 
constant (assuming is classically oscillating, in which case it can be regarded as matter, 
or (f) has decayed into radiation). Relative importance of the effect of cf) compared to that 
of inflaton in ( 1A.26I) depends on the details such as how steep the potential is or how much 
time has passed between horizon exit and the end of inflation. 
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